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Abstract 

We determine the masses and decay constants of pseudoscalar mesons 
D, D s , and K in quenched lattice QCD with exact chiral symmetry. 
For 100 gauge configurations generated with single-plaquette action at 
/3 = 6.1 on the 20 3 x 40 lattice, we compute point-to-point quark propa- 
gators for 30 quark masses in the range 0.03 < m q a < 0.80, and measure 
the time-correlation functions of pseudoscalar and vector mesons. The 
inverse lattice spacing a" 1 is determined with the experimental input of 
f n , while the strange quark bare mass m s a = 0.08, and the charm quark 
bare mass m c a = 0.80 are fixed such that the masses of the correspond- 
ing vector mesons are in good agreement with 0(1020) and J/ip(3097) 
respectively. Our results of pseudoscalar-meson decay constants are 
f K = 152(6) (10) MeV, f D = 235(8) (14) MeV, and f Da = 266(10) (18) 
MeV. 
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1 Introduction 



The pseudoscalar-meson decay constants (e.g., fo, /d s , Ib and fs 3 ) play an 
important role in extracting the CKM matrix elements (e.g., the leptonic decay 
width of — ► l + i>i is proportional to /|> | V^l 2 ), which are crucial for testing 
the flavor sector of the standard model via the unitarity of CKM matrix. 
Experimentally 1 , precise determination of fr> s and fjj will soon result from 
the high-statistics program of CLEO-c, however, the determination of and 
fs s remains beyond the reach of current experiments. 

Theoretically, lattice QCD provides a solid framework to compute the 
masses and decay constants of pseudoscalar mesons (as well as other physical 
observables) nonperturbatively from the first principles of QCD. Thus reliable 
lattice QCD determinations of /# and fg a are of fundamental importance, in 
view of their experimental determinations are still lacking. Obviously, the first 
step for lattice QCD is to check whether the lattice determinations of fo and 
fo s will agree with those coming soon from the high-statistics charm program 
of CLEO-c. This motivates our present study. 

In this paper, we compute quenched quark propagators for 30 quark masses 
in the range 0.03 < m q a < 0.80, in the framework of optimal domain-wall 
fermion proposed by Chiu jT0]-J21- Then we determine the inverse lattice 
spacing a -1 = 2.237(76) GeV from the pion time-correlation function, with 
the experimental input of pion decay constant f n = 132 MeV. The strange 
quark bare mass m s a = 0.08 and the charm quark bare mass m c a = 0.80 
are fixed such that the corresponding masses extracted from the vector meson 
correlation function agree with 0(1020) and J/ip(3097) respectively. Then the 
masses and decay constants of any hadrons containing c, s, and u(d) quarks 2 
are predictions of QCD from the first principles, with the understanding that 
chiral extrapolation to physical m Ut d — m s /25 (or equivalently m n = 135 MeV) 
is required for any observables containing u(d) quarks. 

For pseudoscalar and vector mesons, we measure their time correlation 
functions for the following three categories: (i) two quarks have the same mass; 
(ii) one quark mass is fixed at m s ; (iii) one quark mass is fixed at m c . Note that 
for mesons which are composed of strange and/or charm quarks, their masses 
and decay constants can be measured directly without chiral extrapolation. 

The outline of this paper is as follows. In section 2, we outline our formu- 
lation of exact chiral symmetry on the lattice, and our computation of quark 
propagators. In section 3, we determine the inverse lattice spacing spacing, 
the strange quark bare mass, and the charm quark bare mass. In section 4, we 
present our results of txik and fx- In section 5, we present our results of m D , 
m D s , Id, an d /d s - In section 6, we summarize our results and conclude with 
some remarks. 

1 See, for example, Refs. [Q [21 131 ID > an d other experimental results complied by PDG [5]. 
2 In this paper, we work in the isospin limit m u = m^. 
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2 Lattice quarks with exact chiral symmetry 

To implement exact chiral symmetry on the lattice [HI El E] , we consider 
the optimal domain- wall fermion proposed by Chiu JD]-J2]- The action of 
optimal domain-wall fermion can be written as 

Ns+l 

A F = X^( x > s){(uj s D w (x,x') + 5 XtX ')5 ss i 

s,s'=0 x,x' 

+(u s D w (x, x') - d x , x >)(P+8 s > iS -i + P-5 s i iS+1 )}ij { x', s') (1) 
with boundary conditions 

P + ip(x, —1) = — r m q P + ip(x, N s + 1), 

P-ip(x, N s + 2) = —r m q P^tp(x, 0), r = , 

2m 

where m q is the bare quark mass, and {u s , s = 0, • • • , iV s + 1} are specified 
by the exact formula derived in Ref. [TO] - Here H w = 75-D™, and D w is the 
standard Wilson Dirac operator plus a negative parameter — m (0 < m < 2). 
The quark fields are constructed from the boundary modes at s = and 
s = N s + 1 with ujq = u>n s +i = [12]: 

q( x ) = ^ [P_xp(x, 0) + P + iP(x, N s + 1)] , (2) 
q(x) = V^[^(x,0)P + + ij(x,N s + l)pJ\ . (3) 

After introducing pseudofermions with m q = 2mo, the generating functional 
for n-point Green's function of the quark fields can be derived as [T2*] . 

_ J[dU]e- A °det[(D c + m q ){\ + rD c )- 1 ] exp { J(D C + m q )- l j) 
Z[J ' " = J[dU}e- A sdet[(D c + m q ){l + rD^ 1 ] ^ 

where A 9 is the action of the gauge fields, J and J are the Grassman sources 
of q and q respectively, and 

c - 1 - n 8 =i t s . . 

opt i+nf=iT s ' [) 

rp _ 1 — U) S H W 
1 + O^i^ 

Using the exact formula of uj s JTU], one immediately obtains 

/ H w Rf n \Hl), jV. = 2n + l, 
^ 1 F^r iin) (^), AT s = 2n, W 



where Rz(H^) is the Zolotarev optimal rational polynomial ^3] for the inverse 
square root of H 2 , 



R 



(n,n) / tj2 



do A 1 + 



mm j 

i -r h 2 w /c 2 i-i 



-(/i^ + C 2 „) T~2 | 5 ^tu — #tu/ -^mm (9) 



A m m ™ i = i + C2i_l 



and 



p(n-l,n)/rr2\ _ ^0 n"=l 1 (l + h 2 w /c' 2 i) _ 1 ^ , 1 „s 

A m in {U=l\ L ,L wl c 2l-l) A min t=1 'l> w T <-2«-l 

where the coefficients d , d' , q and are expressed in terms of elliptic functions 
[TB] with arguments depending on N s and A^ oa ./ A^ in , and A mm (A max ) is fixed 
to be the greatest lower bound (least upper bound) of the eigenvalues of \H W \ 
for the set of gauge configurations under investigation. 

From (JU, the effective 4D lattice Dirac operator for the fermion determi- 
nant is 

D{m q ) = {D c + m q ){l + rD.)- 1 = m q + (m - m q /2) [l + l5 H w R z (H 2 w )] (11) 
and the quark propagator in background gauge field is 

5 2 Z[J, J] 



(q(x)q(y)} = - 



j=j=o 



5J(x)5J(y) 

(D c + m q y x l = (1 - rm q y l [D-l{m q ) - r5 x J (12) 



Note that D c is exactly chirally symmetric (i.e. -D c 7 5 + 75-D c = 0) in the 
limit N s — > oo, and its deviation from exact chiral symmetry due to finite 
N s is the minimal provided that the weights {uj s } are fixed according to the 
formula derived in Ref . JU] • Further, the bare quark mass m q (whether heavy 
or light) in the quark propagator (D c + m q )~ l is well-defined for any gauge 
configurations. 

In practice, we have two ways to evaluate the quark propagator (fT2|) in 
background gauge field: 

(i) To solve the linear system of the 5D optimal DWF operator; 

(ii) To solve D~^ y {m q ) from the system 



D{m q )Y 



m q + (mo - m q /2) (l + lh H w R z {H 2 w ))] Y = 1 , (13) 



with nested conjugate gradient [14], and then substitute the solution vector Y 
into 

Since either (i) or (ii) yields exactly the same quark propagator, in princi- 
ple, it does not matter which linear system one actually solves. However, in 
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practice, one should choose the most efficient scheme for one's computational 
system (hardware and software). For our present system (a Linux PC cluster 
of 100 nodes ^S]); it has been geared to the scheme (ii), and it attains the 
maximum efficiency if the inner conjugate gradient loop of (|13|) is iterated with 
Neuberger's 2-pass algorithm [TE]. So we use the scheme (ii) to compute the 
quark propagator, with the quark fields ©-© defined by the boundary modes 
at s = and s = N s + 1. Note that Neuberger's 2-pass algorithm not only 
provides very high precision of chiral symmetry with fixed amount of memory, 
but also is faster than the single pass algorithm for n > 12 ~ 25 (where n is 
the order of the rational polynomial _R( ra_1 > n )) for most computer platforms, as 
discussed by Chiu and Hsieh [T7j . 

We generate 100 gauge configurations with single plaquette gauge action 
at (3 — 6.1 on the 20 3 x 40 lattice. Fixing m = 1.3, we project out 16 low- 
lying eigenmodes of \H W \ and perform the nested conjugate gradient in the 
complement of the vector space spanned by these eigenmodes. For N s = 128, 
the weights {u s } are fixed with X min = 0.18 and \ ma x = 6.3, where X m in < 
^(\H W \) < Xmax for all gauge configurations. For each configuration, point to 
point quark propagators are computed for 30 bare quark masses in the range 
0.03 < m q a < 0.8, with stopping criteria 10~ n and 2 x 10~ 12 for the outer 
and inner conjugate gradient loops respectively. Then the norm of the residual 
vector of each column of the quark propagator is less than 2 x 10~ n 

\\(D c + m q )Y- I|| < 2 x 10~ u , 

and the chiral symmetry breaking due to finite N s is less than 10 -14 , 

Y^S 2 Y 



a 



1 



< lO" 14 , 



for every iteration of the nested conjugate gradient. Further details of our 
scheme have been described in Refs. [TBI ITTj . 

In this paper, we measure the time-correlation functions for pseudoscalar 
(PS) and vector (V) mesons, 

C PS (t) = (£tr{7 5 p c + mQ )-^D c + m q )^ x }\ (14) 

\ x I u 

C v(t) = (^EE tr {7Mp c + ^Q);,o7Mp C + ^)oi}) (15) 

\ 6 /Lt=l x I jj 

where the subscript U denotes averaging over gauge configurations. Here 
Cps(t) and Cy(i) are measured for the following three categories: 

(i) Symmetric masses itlq = m q , 

(ii) Asymmetric masses with fixed itlq = m s = 0.08a _1 , 

(iii) Asymmetric masses with fixed mq = m c = 0.80a _1 , 
where m q is varied for 30 masses in the range 0.03 < m q a < 0.80. 
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3 Determination of a 1 , m s , and m c 
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Figure 1: The pion decay constant f n a versus the bare quark mass m q a. The 
solid line is the linear fit. 



For symmetric masses tuq = m q , the pseudoscalar time-correlation function 
C n (t) ()14|1 is measured, and is fitted to the usual formula 



l^—rrwat _|_ g — m x a(T— t)l 



2m 7r a 

to extract the pion mass m^a and the pion decay constant 

f o ^ 
j^a = 2m q a 



(16) 



1 2 2' 



(17) 



In Figs. H] and 121 we plot the decay constant f n a and pion mass square 
(m^a) 2 versus bare quark mass m g a, respectively. 
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Figure 2: The pion mass square (m^a) 2 versus the bare quark mass m q a. The 
solid line is the fit of Eg. (IT3|) . 

The data of f n a (see Fig. [TJ is well fitted by the straight line 

f n a = 0.059(2) + 0.235(38) x (m q a) . 

Then taking f n a at m q a = equal to 0.132 GeV times the lattice spacing a, 
we can determine the lattice spacing a and its inverse, 

132 

a' 1 = —— GeV = 2.237(76) GeV , 
Jo 

a = 0.088(3) fm . (18) 

Thus the size of our lattice is about (1.8 fm) 3 x 3.6 fm. Since the smallest pion 
mass is 439 MeV, the lattice size is about (3.9) 3 x 7.8, in units of the Compton 
wavelength (~ 0.45 fm) of the smallest pion mass. 

The data of m 2 (see Fig. |2J) can be fitted by the form [T9"j 

m\a 2 = A 1 (m q a)~^ + B{m q a) 2 (19) 
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in quenched chiral perturbation theory (q^PT). The fitted parameters are 



S = 0.187(21) (20) 
Ai = 0.669(45) (21) 
B = 2.666(357) (22) 

with x 2 /d.o.f.=0.54. Evidently, the coefficient of quenched chiral logarithm 
5 = 0.187(21) is in good agreement with the theoretical estimate 5 ~ 0.176 in 
qxPT. 

The bare mass of strange quark is determined by extracting the mass of 
vector meson from the time-correlation function Cy(t). At m q a = 0.08, mya = 
0.460(4), which gives my = 1029(10) MeV, in good agreement with the mass 
of 0(1020). Thus we take the strange quark bare mass to be m s a = 0.08. 
Similarly, at m q a = 0.80, m v a = 1.368(2), which gives m v = 3060(5) MeV, in 
good agreement with the mass of 7/^(3097). Thus, we fix the charm quark 
bare mass to be m c a = 0.80. 



4 fx and rnx 

Next we measure the time-correlation function of kaon Cx(t) (fTf)) with mq 
fixed at m s = 0.08a -1 , while m q is varied for 30 masses in the range 0.03 < 
m q a < 0.80. Then the data of Cx{t) is fitted by the formula analogous to (|T5|l 
to extract the kaon mass mxa and the kaon decay constant /#a. 

In Fig. El the kaon mass nix is plotted versus m„ for 15 quark masses in 
the range 0.03 < m q a < 0.10. The data of m^a can be fitted by 

m K a = 0.197(1) + 0.255(4)(m^a) + 0.389(8) (rn^a) 2 . 

At the physical limit m„ = 135 MeV, it gives rnx = 478(16) MeV, in good 
agreement with the experimental value of kaon mass (495 MeV). 

In Fig. HI fxa is plotted versus bare quark mass m q a. The data is well 
fitted by the straight line 

f K a = 0.068(0) + 0.116(1) x (m q a) 

At m q a = 0, it gives fx = 152(6) MeV, in agreement with the value fx+ = 
159.8 ± 1.4 ± 0.44 MeV complied by PDG 0. 



5 jui fD a , rn D , and m Ds 

Now we turn to charmed pseudoscalar mesons. We measure the time-correlation 
function Cn(t) (J 14)) with mq fixed at m c = 0.80a" 1 , while m q is varied for 30 
different masses in the range 0.03 < m q a < 0.80. Then the data of Cr,(t) is 
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Figure 3: The kaon mass mj{ versus the pion mass m n for 15 bare quark masses 
in the range 0.03 < m q a < 0.10. The solid line is the quadratic fit. 



fitted by the formula analogous to (fTBj) to extract the mass mr>a and decay 
constant f^a. 

In Fig. @ mfla is plotted versus m n a, for 15 quark masses in the range 
0.03 < m q a < 0.10. The data of m^a can be fitted by 

m D a = 0.816(0) + 0.101(3)(m 7r a) + 0.298(6) {m n a) 2 

At m n = 135 MeV, it gives uid = 1842(15) MeV, in good agreement with the 
mass of D meson (1865 MeV). In Fig. the decay constant /do is plotted 
versus bare quark mass m q a. The data is well fitted by the straight line 

f D a = 0.105(1) + 0.172(1) x (m q a) 

At m q a = 0, it gives fo = 235(8) MeV, which serves as a prediction of lattice 
QCD with exact chiral symmetry. 
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Figure 4: The kaon decay constants fxa versus the bare quark mass m q a. The 
solid line is the linear fit. 



The pseudoscalar meson of cs or sc corresponds to mqa = m c a = 0.80 and 
m q a = m s a = 0.08. Its mass and decay constant are extracted directly from 
the time-correlation function, which are plotted as the eleventh data point 
(counting from the smallest one) in Figs. 03 and |Hl respectively. The results 
are m Dg a = 0.878(2) and /u s a = 0.119(2). The mass gives mo s = 1964(5) 
MeV, in good agreement with the mass of _D S (1968). The decay constant gives 
fo s = 266(10) MeV, which agrees with the value f D + = 267±33 MeV complied 
by PDG 0. 



6 Summary and Concluding Remarks 

In this paper, we have determined the masses and decay constants of pseu- 
doscalar mesons K, D and D s , in quenched lattice QCD with exact chiral 
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Figure 5: The mass of D meson m^a versus the pion mass m n a for 15 bare 
quark masses in the range 0.03 < m q a < 0.10. The solid line is the quadratic 
fit. 



symmetry. Our results are: 

478 ± 16 ± 20 MeV, 
1842 ± 15 ±21 MeV, 
1964 ± 5 ± 10 MeV, 
152 ± 6 ± 10 MeV, 
235 ± 8 ± 14 MeV, 
266 ± 10 ± 18 MeV, 

where in each case, the first error is statistical, while the second is our crude 
estimate of combined systematic uncertainty. It is interesting to see whether 
the values of fo and fo 8 coming soon from the high-statistics charm program 
of CLEO-c would agree with our values determined by lattice QCD with ex- 
act chiral symmetry. Further, we note that in a recent 3-flavor unquenched 
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Figure 6: The D-meson decay constants fo-a versus the bare quark mass m q a. 
The solid line is the linear fit. 



lattice QCD study [20] with 0(a 2 ) improved staggered light quarks and 0(a)- 
improved charm quark, their results of fjj and fjj s agree with our values. 

Obviously, our next task is to determinate fs and fs s which are of fun- 
damental importance, in view of their experimental determinations are still 
lacking. Since we will not use any approximations for the heavy b quark, our 
lattice spacing must be small enough such that m^a < 1. Even though this 
does not seem to be formidable for fs s , it is unclear whether we can determine 
Jb reliably via chiral extrapolation. Presumably, /# would behave like a func- 
tion linear in m q for a wide range of m q , similar to fo (FigEJ) and fx (Fig. EJ), 
then one should be able to obtain a reliable chiral extrapolation even for data 
points with m q > m s . 
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